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BETATRON-SYNCHROTRON  DETRAPPING  IN  A 
TAPERED  WIGGLER  FREE  ELECTRON  LASER 


It  has  been  pointed  out  by  M.  Rosenbluth^  that  electron  detrapping  in  the 

ponderomotive  wave  can  occur  if  a  resonance  between  the  electron's  synchrotron 

and  betatron  oscillations  exist.  Synchrotron  oscillations  are  due  to  the 

2-7 

trapped  electrons  oscillating  longitudinally  in  the  ponderomotive  wave, 
while  the  betatron  oscillations  are  predominantly  transverse  electron 
oscillations  due  to  the  transverse  spatial  gradients  associated  with  the 
wiggler  field.®  ^  If  the  radiation  wavefronts  are  curved  these  two 
oscillations  can  be  resonantly  coupled  and  can  lead  to  electron  detrapping  in 
the  ponderomotive  wave.*'^-^®  When  the  radiation  wavefronts  are  curved  the 
electrons  which  are  undergoing  transverse  betatron  oscillations  will 
experience  periodically  different  phases  in  the  ponderomotive  wave.  When  the 
periodically  changing  phase,  due  to  the  transverse  betatron  oscillations,  is 
resonant  with  the  electron's  synchrotron  oscillations  these  later  oscillations 
can  be  amplified  and  result  in  detrapping. 

A  very  qualitative  model  for  this  process  can  be  seen  by  considering  the 

i  A— 7  i  *5-]  ft 

pendulum  equation  •  ’  for  electrons  in  a  wiggler  field  with 

transverse  spatial  gradients  and  a  radiation  field  with  a  curved  wavefront. 

For  electrons  deeply  trapped  in  the  ponderomotive  wave  the  pendulum  equation, 
roughly  speaking,  reduces  to  a  driven  harmonic  oscillator  equation.  The 
characteristic  frequency  of  the  electron  oscillation  is  the  sychrotron 
frequency  and  is  proportional  to  the  fourth  root  of  the  radiation  field 
power.  The  amplitude  of  the  periodic  driving  term  in  the  oscillator  equation 
Is  proportional  to  the  Inverse  of  the  radius  of  curvature  associated  with  the 
radiation  field  and  the  period  of  the  driving  term  is  proportional  to  the 

electron  betatron  period.  The  electron  phase,  described  by  the  pendulum 
Manuacript  approved  October  22,  1984. 
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equation,  can  be  amplified  if  the  frequency  of  the  driving  term  is  resonant 
with  the  characteristic  trapping  frequency.  This  detrapping  mechanism  could 
limit  the  level  of  radiation  power  generated  by  the  FEL. 

In  this  paper  we  suggest  and  analyze  an  alternative  mechanism  for 
sychrotron-betatron  resonant  detrapping  which  does  not  depend  on  the  curvature 
of  the  radiation  wavefronts.  We  show)  that  even  for  a  one  dimensional 
radiation  field  sychrotron-betatron  resonant  detrapping  can  take  place  for 
tapered  wiggler  fields,  i.e.,  the  wiggler  magnetic  field  amplitude. 

To  analyze  this  synchrotron-betatron  resonant  detrapping  mechanism  we 
choose  a  tapered  linearly  polarized  wiggler  field  described  by  the  vector 
potential 


A  (y,z)  =  A  (z)coshfk  (z)ylcosff  k  (z')dz'le  , 
~w  w  w  '  •  w  x 

o 


(1) 


where  A^z)  and  kw(z)  are  the  spatially  slowly  varying  amplitude  and 

wavenumber.  It  will  be  assumed  later  that  kw(z)y  is  somewhat  less  than  unity 
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so  that  cosh(k  y)  «■  1  +  (k  y)  / 2.  The  one  dimensional  radiation  field  is 
w  w' 

described  by  the  vector  potential 


A^(z,t)  =  A^  sin(kz-o)t+0  )e^ , 


(2) 


where  the  amplitude  AL,  wavenumber  k  =  u/c  ,  frequency  w  and  phase  6  are 

assumed  constant.  Since  (1)  and  (2)  are  independent  of  the  x  coordinate,  the 

electron's  canonical  momentum  in  the  x  direction  is  conserved,  i.e., 

d(P  -  lei (A  +  A.)/c)  •  e  /dt  *  0  where  P  *  ym  v  *  (P  ,P  ,P  )  denotes  the 
~  1  1  ~w  ~L  x  ~  o~  x  y  z 

electron's  mechanical  momenta.  The  relativistic  particle  orbit  equations 
become 
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where  A(y,z,t)  =  fA^Cy.z)  +  A^(z,t)l*  ex  is  the  total  x  component  of  the 
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vector  potential  and  y  *  (1  +  P  •  P/mQc  ) *  •  In  obtaining  (3)  we  used  the 

fact  that 
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and  assumed  chat  the  injected  momentum  in  the  x  direction  is  zero,  i.e., 
Pqx  =  0.  Using  (3)  and  (4)  the  electron's  axial  velocity,  v^,  is  given  by 
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A  (y,z,t). 
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In  what  follows  it  proves  convenient  to  perform  a  transformation  from  the 
independent  time  variable  t  to  the  independent  position  variable  z.  An 
electron's  phase  with  respect  to  the  portderomotive  wave,  in  terms  of  the 
position  variable  z,  is  defined  as 


~  z 

*<VV*o*8>  °  ^o  +  ^  +  k  ~  ta)/vz^y0>'i;0*'i,o»z'^)dz'»  (ft) 

o 


where  yQ  is  the  electron's  y  position  at  z  "  0,  i|<0  is  the  phase  at  z  *  0, 

ij»'  *  3 <p/ 3z  I  =  k  (0)  +  k  -  u/v  and  v  is  the  axial  velocity  at  z  ■  0 
o  _  w  o  o 

z30 

(assumed  to  be  identical  for  all  electrons).  Quantities  denoted  with  the 


superscript  ~  are  functions  of  the  initial  condition  variables  and 
the  independent  variable  z.  Thus,  for  example,  v^  is  the  axial  velocity  at 
position  z  of  an  electron  with  initial  condition  (at  z  *  0)  variables 
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yo,^o  and  Differentiating  (6)  twice  yields 
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where  '  denotes  the  operator  3/3z.  Upon  performing  the  operation 
2 

3/ 3z  +  v^c  3/ 3t  in  (5),  transforming  the  resulting  equation  from  the 
independent  variable  t  to  the  variable  z  and  substituting  the  result  into  (7), 
we  arrive  at  the  generalized  pendulum  equation  for  a  tapered  wiggler  with 
transverse  spatial  gradients, 
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r  f  2A  k'  cosh  (k  y)  +  A  k'y  sinh(2k  y) )cos  ff  k  dz') 
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-  k  A  cosh  (k  y)  sin  (2  f  k  dz") 
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+  ( A'k  cosh(k  y)  +  A  ATk"y  sinh(k  y)  )sim{/ 
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(k^  +  k  -  v^w/c  )A^A^ cosh (k^y) cos tyl 


(8) 


where  v^  =  “/(k^  +  k  -  if)')  and  y  =  y(z).  The  last  term  in  (8)  represents  the 
usual  ponderomotive  potential  wave. 

Before  going  on  to  simplify  (8)  an  equation  describing  the  electrons 
betatron  motion  (transverse  oscillations  in  y)  is  needed.  Using  (3a)  and 
assuming  | k^y | <<1  and  | A^ | << | A^ |  we  find  that 
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In  arriving  at  (9)  we  have  also  assumed  that  |yy|<<yy,  these  approximations 
can  be  shown  to  be  well  satisfied.  Transforming  (9)  from  the  independent 
variable  t  to  z  and  setting  v^  =*  c  yields  the  following  equation  for  the 
betatron  orbits. 


y"  +  k  (z)fl  +  cos(2  f  k  (z')dz"))y  *  0, 

p  w 


(10) 


a 

where  k„(z)  =  1  e  I A  k  /(/  2  yin  c  )  =*  8  k  / /2  and  6  ■  I  e  I A  (z)/ym  c  is  the 
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normalized  wiggle  velocity.  Since  the  betatron  wavenumber,  k  ,  is  much 
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greater  than  kw,  the  betatron  motion  can  be  separated  into  a  slowly  varying 
part  and  a  small  rapidly  varying  part.  Neglecting  the  rapid  variations 
in  y(z),  the  solution  of  (10)  is 


y(z)  =  y(0)  (k  (0)/k  (z)  ^2cos  (  f  k  (z')dz'  +  4>(0)),  (11) 

P  P  P 

o 


where  y(0)  and  <J>(0)  are  constants. 

We  can  now  proceed  to  simplify  the  pendulum  equation  in  (8).  Assuming 
~  2  2 

(k  y)  «  1,  A'  «  k  A  ,  k'  «  k  and  keeping  only  the  slowly  varying  terms, 
w  w  w  w  w  w 

(8)  reduces  to  the  form 
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where  Kg(z)  “  (|e|  wk^A^A^/fY  mQc  ))  Is  the  synchrotron  wavenumber, 

a2  -  | e | 2w/(4Y2m^c^) ,  £(z)  •  y2(0)  (kfl(0)/2ko(z) )d(A  k  )2/dz  and  (11)  was  used 

0  P  p  w  w 

to  replace  y(z).  We  can  further  simplify  (12)  by  noting  that 
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Therefore  we  may  keep  only  the  driving  term  associated  with  the  betatron 

oscillations  which  could  amplify  the  synchrotron  oscillations,  l.e., 
z 

cos (2  /  k  dz'  +  2<p(0)],  and  (12)  reduces  to 
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K^(z)sint{/  = 
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dk 
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-  a  +  C(z)cos(2  f  kRdz  ')  ] , 
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(13) 


where  for  later  convenience  we  have  shifted  the  phase  ip  by  n/2, 

l.e.  ip  =  ip  +  n/2,  and  set  <{>(0 )  =  0.  Note  that  the  electron's  energy, 

neglecting  transverse  wiggler  gradients,  is  determined  by  the  relation 

I e I A  k 

9y/3z  - - -  ■  yffi-  ■  sinip.  (14) 

2Ymoc 

Defining  the  resonant  phase,  ipR,  in  the  usual  way,  l.e.. 


sinipR 


fdkw/dz  -  a2dA2/dz)/K2, 


(15) 


the  pendulum  equation  (13)  becomes 
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— f  +  K  sinijj  =  K  (l  +  e  cos(2  f  k adz ')  )sini|L 
dz  o 

82  dk  z 

-  e(l  +  2^  £“ )  cos  (2/  kgdz'), 

w  o 

where  e  =  k2  y2(0)/2. 
w 

We  now  consider  the  particularly  simple  illustration  of  a  wiggler  field 

with  a  constant  period  Ok  /3z  =  0)  and  linearly  changing  amplitude, 

w 

A  (z)  =  A  (0)  +  6A  z/L  ,  where  A  (0)  »  I 6A  I  is  constant  and  L  is  the  length  of 
w  w  w  w  w  1  w1  w 

the  wiggler  field.  The  pendulum  equation  in  (13)  for  the  case  where 

|6Aw/Aw(0)|  «  1,  Kg  =  2kg/(l  +  6)  and  |s|  «  1,  i.e.,  the  synchrotron  and 

betatron  oscillations  are  resonant,  reduces  to 
.2 

— ~  +  sintfi  =  (l  +  e  cos(l  +  6)z)sintJ>  ,  (17) 

dZ2  R 

2 

where  Z  =  K  z  and  sinit  =  -  y  f 6A  /A  (0))/(2k  L  ).  For  the  case  where  A,,  is 
s  Rz  ww'ww  w 

constant  and  k^(z)  =  kw(0)  +  6kwz/Lw  var*es  linearly,  the  pendulum  equation  in 
(13)  reduces  to 


+  sini(i  =  fl  -  tP2y2  cos ( 1  +  6)z)siniJ>  ,  (18) 

dZ  w  z 

where  sinil)  =  (6k  /k  (0)/(282L  k  (0)). 
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Rather  than  obtain  an  approximate  solution  to  (17)  or  (18),  using  a  multiple 
time  scale  approach,  we  simply  solve  (17)  numerically.  Initially  the  phases  are 
distributed  uniformly  between  0  and  2tt  with  3ip/ 3z  =  0.  Figure  (1)  shows  the 
precentage  of  trapped  particles  as  a  function  of  normalized  distance  for 
sinijL  =0.3  and  e  =  0. 1  and  0.15.  Initially  approximately  65%  of  the  particles 
are  trapped  and  for  e  =  0  this  fraction  is  of  cause  maintained.  After  about  5 
synchrotron  oscillations  55%  are  trapped  for  e  =  0.1  and  50%  for  e  =  0.15. 


Since  we  have  assumed  zero  beam  emittance,  these  results  apply  only  to  those 


electrons  Initially  on  the  outer  edge  of  the  beam,  i.e.,  those  having  the  largest 
value  of  e.  Figure  (2)  shows  the  percentage  of  particles  detrapped  after  10 
synchrotron  oscillations  as  a  function  of  the  mismatch  parameter  6,  for 
sin<|/_  =  0.3  and  e  =0.1,  0.15.  The  percentage  of  detrapped  particles  maximize 
when  6  <  0  since  the  more  deeply  trapped  particles  oscillate  slightly  faster  than 
those  trapped  nearer  the  phase  space  sepratrlx.  Here  again  these  results  apply 
to  only  a  small  fraction  of  the  total  number  of  beam  electrons,  those  initially 
near  the  edge  of  the  beam. 

Our  model  is  somewhat  Idealized  since,  among  other  things,  beam  emlttance 
has  been  neglected.  The  neglect  of  emlttance  implies  that  all  the  electrons  have 
the  same  Initial  betatron  phase,  $(0),  see  Eq.  (11).  Hence,  those  electrons 
Injected  near  the  axis  will  not  experience  the  betatron  synchrotron  detrapping 
since  the  value  of  e  for  these  electrons  Is  much  smaller  than  those  near  the  edge 
of  the  beam. 
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Fig.  1  —  Percentage  of  trapped  electrons  as  a  function  of  normalized 
distance  for  sini/zj^  =  0.3,  and  two  values  of  e. 
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Fig.  2  —  Percentage  of  trapped  electrons  after  10  synchrotron  oscillations  as  a  function  of 
the  mismatch  parameter  5  for  sin^R  =  0.3  and  e  =  0.1,  and  0.15. 
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